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Abstract 

A  novel  radial  flow  packed-bed  microreactor  has  been  developed  to  enable  hydrogen  delivery  in  miniature  fuel  cells  for  portable  power 
generation  using  the  methanol  reforming  process.  Conventional  rectangular  microchannel  based  designs  typically  used  for  on-chip  chemical 
processing  were  found  to  be  inadequate  for  practical  high-throughput  micro  packed-bed  application  due  to  the  associated  high  pressure  drops. 
The  radial  design  exploits  several  physical  phenomena  inherent  to  the  device  operation  and  provides  an  order  of  magnitude  reduction  in  the 
required  driving  pressure  while  maintaining  the  compact  design  of  the  microreactor  and  providing  greater  than  98%  conversion  of  methanol 
to  hydrogen  at  flow  rates  sufficient  to  supply  hydrogen  to  a  20  W  fuel  cell.  The  developed  microreactor  opens  the  path  for  fuel  cell  based 
portable  energy  systems  with  energy  storage  capacities  up  to  10  times  those  possible  with  current  battery  technology. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction  and  literature  survey 

Miniaturized  Tab-on-a-chip’  and  ‘plant-on-a-chip’  de¬ 
vices  have  been  successfully  developed  for  several  real-world 
applications  in  recent  years  [1-4].  Power-plant-on-chip  and 
palm  power  are  buzzwords  of  great  interest  currently  in  the 
micro-power  community  where  battery  technology  is  thought 
to  be  quickly  approaching  its  peak  in  energy  storage  capacity 
per  unit  volume/weight  of  the  power  source  [5].  As  a  promis¬ 
ing  alternative,  hydrogen  fuel  cells  have  long  been  considered 
a  power  source  of  the  future  due  to  their  high  efficiency  in  the 
conversion  of  fuel  to  usable  energy  and  low  emission  of  pollu¬ 
tants  and  greenhouse  gases  during  operation  [6,7].  In  the  con¬ 
text  of  portable  and  palm  power  sources,  fuel  cells  have  the 
potential  of  providing  energy  storage  densities  several  times 
those  possible  using  current  state-of-the-art  lithium-ion  bat¬ 
teries  [5].  However,  the  difficulties  and  hazards  involved  in 
the  storage  and  handling  of  the  hydrogen  fuel  in  either  com¬ 
pressed  gas  or  liquid  form  have  been  a  major  impediment  in 
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the  successful  commercialization  of  miniature  hydrogen  pro¬ 
ton  exchange  membrane  (PEM)  fuel  cells.  Also,  the  stored 
energy  density  of  hydrogen  in  compressed  or  liquid  form  is 
significantly  lower  in  comparison  to  storage  in  the  form  of  liq¬ 
uid  hydrocarbons,  such  as  methanol,  which  can  be  reformed 
for  generation  of  the  hydrogen  gas  as  and  when  needed.  Hy¬ 
drogen  storage  in  metal  hydrides  and  nano  structured  carbon 
[8,9]  has  been  proposed  but  a  number  of  limitations  have  been 
recognized,  including  high  weight  of  the  hydride  storage  unit, 
lower  energy  density,  the  passive  nature  of  the  hydrogen  re¬ 
lease  mechanism  and  high  cost.  Therefore,  the  development 
of  a  compact  and  efficient  miniaturized  fuel  ‘reformer’  for 
on-demand  hydrogen  delivery  to  miniature  PEM  fuel  cells 
has  been  the  subject  of  considerable  research  activity  in  re¬ 
cent  years  [10-19]. 

2.  The  radial  microreactor 

In  this  paper,  we  present  a  novel  radial  flow  micro  packed 
bed  reactor-on-a-chip  with  an  integrated  micro-vaporizer  for 
hydrogen  delivery  using  methanol  reforming: 

CH3OH  +  H20  ^  C02  +  3H2  (1) 
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CH3OH  ^  CO  +  2H2  (2) 

C0  +  H20^C02  +  H2.  (3) 

We  have  previously  demonstrated  rectangular  microchannel 
based  packed-bed  microreactors  for  this  application  [19]  as 
part  of  our  ongoing  work  on  this  subject.  The  packing  of  mi¬ 
crochannels  (400  |jim  x  1000  fxm)  with  catalyst  particles  of 
mean  diameter  50-70  p,m  resulted  in  highly  efficient  contact¬ 
ing  of  reactants  with  the  catalyst  and  excellent  conversions 
at  relatively  low  temperatures  in  the  range  of  195-210  °C. 
However,  the  inlet  pressure  of  the  methanol-water  mixture 
required  to  drive  the  flow  and  generate  sufficient  hydrogen 
for  an  8-10  W  fuel  cell  was  routinely  found  to  be  in  the  range 
of  70-100  psig  in  these  conventional  microchannel  reactors. 
Such  high  inlet  pressures  are  neither  achievable  nor  desir¬ 
able  in  portable  devices  nor  can  they  be  controlled  using 
available  miniature  pumping  devices  and  valves  [20,21].  For 
true  high  energy  density  portable  operation,  the  microreactor 
design  should,  therefore,  be  able  to  operate  at  significantly 
lower  pressure  drops  while  maintaining  the  high  throughput 
rate  for  hydrogen  production.  The  unique  radial  design  pre¬ 
sented  here  achieves  this  goal  by  exploiting  several  key  char¬ 
acteristics  of  the  methanol  reforming  reaction  as  discussed 
below. 

An  analysis  of  the  different  physical  and  chemical  phe¬ 
nomena  occurring  inside  the  microchannel  reactor  reveals  the 
inadequacy  of  simple  microchannel  based  designs  for  practi¬ 
cal  micro-packed  bed  applications.  The  methanol  reforming 
reaction  causes  an  increase  in  the  total  number  of  moles  of 
the  reaction  mixture,  leading  to  an  increase  in  gas  mixture 
volume  at  a  given  pressure.  This  results  in  an  increase  in 
the  gas  flow  velocity  towards  the  outlet  of  the  microchan¬ 
nels  where  a  significant  amount  of  methanol  has  already 


been  converted  to  hydrogen.  At  the  same  time,  the  pressure 
decreases  as  the  gas  moves  closer  to  the  outlet  due  to  the 
presence  of  a  pressure  gradient  which  drives  the  flow  and 
this  causes  a  further  increase  in  the  gas  volumetric  flow  rate 
and  gas  flow  velocity  [19].  The  result  of  both  these  effects 
is  that  higher  velocity  in  local  regions  closer  to  the  outlet 
leads  to  a  smaller  residence  time  in  these  regions,  result¬ 
ing  in  lower  catalyst  utilization.  Also,  since  the  bulk  of  the 
pressure  drop  in  the  microchannels  occurs  in  the  high  gas 
velocity  regions,  such  a  design  leads  to  high  pressure  drops 
in  the  device  during  operation,  as  was  experimentally  ob¬ 
served. 

Vaporization  of  liquid  feed  in  a  microreactor  also  poses 
unique  challenges  due  to  the  sudden  expansion  of  the  fluid 
during  boiling  which  causes  flow  pulsation  and  severe  pres¬ 
sure  drop  oscillation  [22]  in  the  two-phase  flow  region.  This 
pulsation  in  the  vaporization  region  leads  to  unsteady  opera¬ 
tion  and  the  possibility  of  damage  to  the  device  unless  suffi¬ 
cient  room  is  provided  for  the  nearly  thousandfold  expansion 
of  the  fluid  as  it  changes  from  liquid  to  vapor. 

All  of  the  above  problems  have  been  resolved  by  employ¬ 
ing  a  radial  flow  microreactor  geometry  which  compensates 
for  the  expected  gas  expansion  by  providing  a  variable  (in¬ 
creasing)  flow  cross  section  along  the  reaction  path.  As  shown 
in  Fig.  1(a),  the  reacting  gases  flow  radially  outward  from  the 
smaller  inner  circular  boundary  to  the  larger  outer  circular 
boundary,  resulting  in  significantly  reduced  pressure  drops 
compared  to  conventional  microchannel  based  designs  for  a 
given  feed  flow  rate.  While  radial  flow  reactors  have  been 
used  in  conventional  chemical  plants  in  the  past  [23],  their 
use  has  been  limited  due  to  the  more  complex  manufacturing 
process  involved  in  implementing  these  designs  at  the  macro 
scale.  This  is  not  an  issue  in  micro  scale  radial  reactors,  since 
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Fig.  1.  Radial  flow  microreactor:  (a)  schematic;  (b)  fabricated  device;  (c)  vaporization  region;  and  (d)  catalyst  particle  filter  (patent  pending). 
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the  microfabrication  steps  involved  are  not  much  different  for 
the  radial  and  microchannel  configurations. 

Apart  from  the  lower  pressure  drops,  the  radial  design 
results  in  a  very  compact  device,  since  a  significantly  larger 
amount  of  catalyst  can  be  accommodated  in  a  reactor  of  given 
volume  compared  to  a  microchannel  based  design  in  which 
the  reactor  walls  take  up  a  significant  portion  of  the  net  reactor 
volume  [19] .  While  this  may  not  always  be  an  important  issue 
in  conventional  chemical  plants  which  are  typically  the  size  of 
large  buildings,  compactness  of  design  is  of  prime  importance 
in  portable  devices  where  any  increase  in  device  size  may 
severely  limit  its  applicability.  Thus,  the  radial  microreactor 
exploits  the  unique  advantages  of  the  radial  design  at  micro¬ 
scale,  making  the  device  especially  well  suited  for  plant-on-a- 
chip  applications  where  low  pressure  drop  and  compactness 
of  design  are  of  prime  importance. 


3.  Fabrication  procedure 
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A  schematic  of  the  steps  followed  in  fabricating  the  ra¬ 
dial  microreactor  are  given  in  Fig.  2.  Standard  photolithog¬ 
raphy  steps  followed  by  deep-reactive  ion  etching  (DRIE) 
are  used  for  etching  features  on  a  silicon  substrate  to  form 
the  vaporizer  and  filters  for  trapping  catalyst  particles.  This 
is  followed  by  anodic  bonding  to  a  pyrex  wafer  having 
drilled  holes  for  catalyst  introduction  (two  catalyst  intro¬ 
duction  ports)  and  microfluidic  interfacing  (one  inlet  and 
two  symmetrically  located  outlets  as  shown  in  Fig.  1(b)). 
The  procedure  for  carrying  out  the  microfluidic  interfacing 
is  described  elsewhere  [24].  Finally,  the  catalyst  is  intro¬ 
duced  into  the  annular  reaction  chamber  by  carrying  out 
fluidized  packing  of  Cu/ZnO/Al2C>3  microparticles  ( Sud - 
Chemie  catalyst  Cl 8-7)  of  diameters  in  the  range  of  50- 
70  |xm.  The  catalyst  powder  is  placed  on  the  catalyst  intro¬ 
duction  ports  and  the  inlet/outlet  ports  of  the  microreactor 
are  connected  to  vacuum.  This  draws  the  microparticles  into 
the  reaction  chamber  where  they  are  trapped  in  the  annular 
region  between  the  catalyst  particle  filters  to  form  the  radial 
packed  bed  shown  in  Fig.  2(e).  The  fabricated  device  along 
with  magnified  views  of  the  different  features  is  shown  in 
Fig.  l(b)-(d). 


4.  Experimental  verification 

4.1.  Experimental  set-up 

A  schematic  of  the  experimental  set-up  used  for  testing  the 
microreactor  is  shown  in  Fig.  3  and  the  actual  set-up  is  shown 
in  Fig.  4.  It  consists  of  a  syringe  pump  for  feeding  the  liq¬ 
uid  methanol-water  mixture  of  required  composition  to  the 
microreactor  at  desired  flow  rates  and  a  mass  spectrometer 
(Stanford  Research  Systems  QMS-200  gas  analyzer)  which 
provides  an  on-line  analysis  of  the  composition  of  the  re¬ 
acted  gases  from  the  microreactor  exit.  An  inert  gas  (argon) 


(c)  Microfluidic  device  after  anodic  bonding 
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(e)  Formation  of  radial  micro  packed  bed 


Fig.  2.  Schematic  of  the  steps  followed  in  the  fabrication  of  the  radial  mi¬ 
croreactor. 


of  known  flow  rate  maintained  using  a  mass  flow  controller 
is  mixed  with  the  microreactor  exhaust  gases  before  feeding 
into  the  gas  analyzer.  This  provides  for  calculation  of  the  ex¬ 
act  hydrogen  production  rate  by  measuring  the  concentration 
of  hydrogen  relative  to  argon.  A  pressure  sensor  connected 
between  the  syringe  pump  and  the  microreactor  inlet  is  used 
to  monitor  the  pressure  of  the  methanol-water  mixture  enter¬ 
ing  the  microreactor. 

The  mass  spectrometer  takes  in  the  desired  amount  of  gas 
for  analysis  and  the  rest  is  either  vented  through  a  vacuum 
pump  connected  to  the  outlet  side  (Fig.  3)  or  collected  in  a 
sample  balloon  as  shown  in  Fig.  4.  The  temperature  of  the 
microreactor  is  maintained  at  desired  levels  using  a  hotplate 
on  which  the  microreactor  is  placed.  As  demonstrated  previ- 
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Fig.  3.  Schematic  of  experimental  setup. 


ously  [19],  it  is  also  possible  to  use  on-chip  heaters  and  tem¬ 
perature  sensors  for  maintaining  the  operating  temperature 
of  the  microreactor  (Fig.  3)  to  avoid  the  use  of  bulky  external 
hotplates  and  simplify  the  experimental  set-up. 

4.2.  Pressure  drop  experiments 

It  was  experimentally  verified  that  this  device  could  be  op¬ 
erated  at  less  than  ^th  the  pressure  drops  and  with  through¬ 
puts  up  to  three  times  those  possible  using  a  microchannel 
configuration  of  the  same  overall  reactor  size,  as  shown  in 
Fig.  5. 

4.3.  Vaporization  experiments 

The  radial  geometry  of  the  design  provides  for  expan¬ 
sion  of  the  evaporating  liquid  radially  outward  from  the  inlet 
at  the  center  and  successfully  dampens  the  pulsation  effect 
discussed  earlier.  Special  features  were  incorporated  in  the 
device  as  shown  in  Fig.  1(c)  to  provide  sufficient  area  of  con¬ 
tact  between  the  liquid  feed  and  the  heated  silicon  surface 
to  ensure  proper  vaporization.  The  central  radial  micro- walls 
in  Fig.  1(c)  serve  as  the  initial  point  of  contact  for  the  liquid 
feed  entering  at  the  center.  Subsequently,  the  two-phase  fluid 
mixture  meanders  around  the  concentric  walls  and  passes  on 
to  the  catalyst  packed  reaction  zone  through  another  set  of 
parallel  walls  providing  further  area  for  vaporization  of  the 


Fig.  5.  Comparison  of  pressure  drops  in  the  microchannel  device  with  pres¬ 
sure  drops  in  the  radial  device  at  three  times  the  throughput  rate  (feed:  2:1 
molar  H20:CH30H). 

liquid  feed.  Though  the  presence  of  the  meandering  circular 
walls  causes  the  feed  to  move  in  non-radial  directions  within 
the  vaporization  region,  the  features  are  arranged  such  that 
the  fluid  exits  the  vaporization  region  uniformly  in  all  di¬ 
rections.  This  was  verified  by  feeding  liquid  water  at  room 
temperature  at  the  inlet  and  noting  that  the  flow  out  of  the 
vaporization  zone  was  uniform  in  all  directions  [25]. 

The  feed  vaporization  ability  of  the  device  was  tested  sepa¬ 
rately  prior  to  carrying  out  the  reaction  runs  by  feeding  liquid 
water  at  various  flow  rates  and  temperatures  [26-28]  and  not¬ 
ing  the  maximum  flow  rate  up  to  which  the  entire  feed  could 
be  vaporized  within  the  vaporization  region.  The  device  was 
able  to  successfully  vaporize  up  to  15-25  cm3  h_1  of  liquid 
water  feed  at  operating  temperatures  in  the  range  of  150- 
200  °C  without  any  liquid  droplets  appearing  in  the  catalyst 
region,  as  shown  in  the  vaporization  runs.  Snapshots  from 
one  of  the  vaporization  runs  at  200  °C  [28]  and  water  feed 
flow  rates  of  15  and  20  cm3  h_1  are  shown  in  Fig.  6(a)  and 
(b),  respectively.  Liquid  droplets  can  be  seen  clearly  within 
the  vaporization  region  at  the  high  20  cm3  h-1  flow  rate,  but 
all  fluid  entering  the  catalyst  zone  (marked  by  the  outer  cir¬ 
cular  boundary  of  the  vaporizer  shown  in  Fig.  6(a)  and  (b)) 
is  in  gaseous  phase  (note  the  absence  of  liquid  droplets  in 
this  region),  thus,  verifying  complete  vaporization.  Since  the 
maximum  flow  rate  of  the  methanol-water  liquid  feed  during 
the  subsequent  reaction  runs  needed  to  produce  hydrogen  for 
a  20- W  fuel  cell  was  less  than  10  cm3  h-1,  the  vaporization 
capacity  was  demonstrated  to  be  sufficient  for  the  range  of 
interest  of  the  feed  flow  rates. 


5.  Results  from  reaction  runs 

Results  from  experimental  reaction  runs  of  the  radial  re¬ 
actor  are  shown  in  Table  1  and  in  Fig.  7.  The  temperature 
of  the  microreactor  was  varied  in  the  range  of  150  to  about 
250  °C  for  different  feed  flow  rates  and  feed  compositions. 
The  hydrogen  production  rate  for  each  run  rises  as  increas¬ 
ing  the  reactor  operating  temperature  results  in  faster  kinetics 
to  the  point  when  equilibrium  conversions  are  reached  and 
no  further  increase  in  hydrogen  production  rate  is  possible. 
Since  the  methanol  reforming  reaction  can  proceed  to  almost 
complete  conversion  at  equilibrium  under  the  operating  con- 


Fig.  4.  The  experimental  setup. 
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ditions  studied  here  [29],  this  represents  close  to  complete 
conversion  of  methanol  in  the  radial  reactor.  The  conversion 
values  were  calculated  using  the  data  reported  in  the  table  and 
were  found  to  be  higher  than  98%  for  each  of  the  reaction 
runs.  As  shown  in  Table  1,  higher  water  to  methanol  ratios 
in  the  feed  lead  to  an  increase  in  the  number  of  moles  of  hy¬ 
drogen  produced  per  mole  of  methanol  in  the  feed  due  to  the 
suppression  of  carbon  monoxide  production  via  reaction  (3) 
above. 

The  peak  hydrogen  production  rate  achieved  with  the  ra¬ 
dial  device  was  about  145  SCCM  (sufficient  for  a  20  W  fuel 
cell),  as  shown  in  Fig.  7.  The  pressure  drop  across  the  de¬ 
vice  at  this  peak  production  rate  was  found  to  be  10.2  psi. 
In  contrast,  the  microchannel  device  developed  earlier  [19] 
provided  a  peak  hydrogen  production  rate  of  72  SCCM,  at  a 
pressure  drop  exceeding  95  psi.  Thus,  it  was  experimentally 
verified  that  the  radial  device  can  produce  more  than  twice  the 
amount  of  hydrogen  compared  to  the  microchannel  device, 
at  less  than  gth  of  the  pressure  drop. 


♦  Feed:  7.0  cc/hr  liquid,  1. 0:1.0  molar  water: methanol 
■  Feed:  7.5  cc/hr  liquid,  1 .5:1 .0  molar  water: methanol 
▲  Feed:  8.0  cc/hr  liquid,  2.0:1 .0  molar  water: methanol 
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Fig.  7.  Hydrogen  production  rate  in  radial  microreactor  at  various  tempera¬ 
tures  and  feeds. 


6.  Mathematical  analysis 

Rigorous  mathematical  modeling  and  analysis  was  carried 
out  to  understand  the  observed  improvement  in  performance 
with  the  radial  configuration  in  comparison  to  the  microchan¬ 
nel  configurations  developed  earlier.  The  main  equations  of 
the  model  are  discussed  below,  followed  by  results  from  nu¬ 
merical  integration  of  the  developed  models. 


6.1.  Coupled  flow  and  reaction  kinetics  modeling 


A  comprehensive  model  based  on  coupled  reaction  kinet¬ 
ics  and  momentum  transport  equations  is  provided  below. 
FORTRAN  programs  were  developed  to  integrate  the  model 
equations  numerically  using  the  method  of  lines.  As  dis¬ 
cussed  below,  the  simulation  results  support  the  experimen¬ 
tally  observed  behavior  (significantly  lower  pressure  drops 
at  higher  throughputs  in  the  radial  device  compared  to  the 
microchannel  device)  and  serve  as  further  verification  of  the 
advantage  of  the  radial  geometry  for  high  throughput,  low 
pressure  drop  applications. 

Darcy’s  law  for  flow  through  porous  media  (Eq.  (4))  [30] 
and  the  continuity  equation  for  compressible  fluids  (Eq.  (5)) 
[3 1  ]  are  used  to  model  gas  flow  through  the  micro-packed  bed, 
where  u  is  the  fluid  flow  velocity  vector;  /Tp,  the  permeability 
of  the  porous  medium;  P  denotes  the  pressure;  /x  denotes  the 
fluid  viscosity  and  p,  the  fluid  density: 

Kv 

u  =  uxi  +  uyj  +  uzk  = - -VP  (4) 

/a 

V  •  (p  x  u)  =  0  (5) 


The  gas  is  modeled  as  a  compressible  fluid  in  accordance 
with  the  ideal  gas  law: 

PV  =  nRT  (6) 


Table  1 


Results  from  experimental  operation  of  the  radial  microreactor 


Feed  composition 
(molar  H20:CH30H) 

Feed  flow  rate 
(liquid,  20 °C)  (cm3!!-1) 

Moles  of  H2  produced 
per  mole  CH3OH  feed 

Mol%  CO 
in  product 

1:1 

7.0 

2.83 

3.13 

1.5:1 

7.5 

2.85 

2.52 

2:1 

8.0 

2.96 

2.18 
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where  V  denotes  the  volume  occupied  by  n  moles  of  the  gas 
at  temperature  T  and  R  is  the  universal  gas  constant.  This  can 
also  be  written  as 

nM  MP 

- =  p  = -  (7) 

V  RT 

where  M  denotes  the  gas  molecular  weight.  Eqs.  (4),  (5) 
and  (7)  can  be  combined  into  a  single  governing  equation 
describing  the  flow  of  the  gas  mixture  through  the  packed 
bed  of  catalyst  particles: 


denoting  the  forward  rate  of  the  reactions  (molm-3s-1) 
given  in  Eqs.  (15)— (17)  above  by  r\,  r2,  and  r3,  respectively 
the  rate  of  production  of  each  species  is  given  by: 


^ch3oh  =  — n  -  r2  (18) 

m2o  =  ~n  -  r3  (19) 

rco2  =ri+r3  (20) 

rn2  =  3ri  +  2r2  +  r3  (21) 


KVM 

PWP 

piRT 

and  since  Kv  and  R  are  constants  and  assuming  that  T  is 
maintained  constant  across  the  reactor  geometry  by  proper 
temperature  control,  this  becomes: 


V  • 


PVPj 


rco  =n~r3  (22) 

The  rates  r\ ,  V2  and  r 3  are  functions  of  the  partial  pressures 
(Pj,  N  m-2)  of  the  different  species  in  the  reacting  mixture, 
which  are  given  by: 

Pj  =  nj  x  P  (23) 

where  P  is  the  total  pressure  (N  m-2)  and  nj  is  the  mole 
fraction  of  each  species,  given  by: 


Considering  variations  only  in  radial  direction  (axisym- 
metry),  we  have: 


(24) 


1  d  f  MPdP\ 

-  —  r - =0 

r  dr  \  pi  dr  J 


mp  dP 


pi  dr 


(10) 

(11) 


rMP  d2P  rM  fdP\2  rP  dP  dM  MP  dP 

yu  dr2  y,  \  dr  J  yu  dr  dr  yu  dr 

rMP  dP  dpi 

- 0 - -  =0. 

piz  dr  dr 

Then  using  the  notation 


(12) 


( 


dP 

dr 


) 


=  m 


(13) 


and  removing  the  partial  derivative  sign  for  this  one¬ 
dimensional  (radial)  representation,  we  have  from  Eqs.  (12) 
and  (13)  above, 

dm 
dr 

Eqs.  (13)  and  (14)  constitute  the  fluid  mechanics  model  for 
the  radial  reactor.  In  this  model,  the  effect  of  reaction  is  cap¬ 
tured  through  variations  in  effective  molecular  weight  ‘AT 
and  viscosity  ‘/x’  of  the  reacting  mixture.  For  the  methanol 
reforming  system, 


nr 


m  m  dM  m  dpi 

—  — - _l_ - 

r  M  dr  pi  dr 


(14) 


CH3OH  +  H20  ^  C02  +  3H2  (15) 

CH3OH^CO  +  2H2  (16) 

CO  +  H20  ^  C02  +  H2  (17) 


where  represents  any  of  CH3OH,  H20,  C02,  CO  and  H2 
and  F  is  the  flow  rate  in  mols-1.  This  is  in  turn  used  to 
calculate  the  molecular  weight  M  for  use  in  Eq.  14  via  the 
equation: 

n 

M  =  njMj.  (25) 

j=  1 


Also,  the  effect  of  variation  in  viscosity  of  the  reaction 
mixture  as  its  composition  changes  with  reaction  is  included 
in  Eq.  (14)  above.  The  mixture  viscosity  (pi),  which  is  a  func¬ 
tion  of  the  temperature  (T),  pressure  (P),  and  composition 
(rij),  is  calculated  during  the  numerical  solution  by  using 
correlations  and  mixing  rules  available  in  the  literature  [32]. 

In  order  to  incorporate  the  effect  of  chemical  reactions 
in  the  above  fluid  dynamics  model,  we  write  the  following 
balance  for  the  flow  rate  Fj  (mol  s_1)  of  each  species  in  the 
radial  flow  device: 


d  F 


1 


d  F 


1 


dV  2nrh  dr 


_ Q  _  1 

=  rj  (mol  (m  s  )) 


(26) 


=  ( 2nrh)rj 


(27) 


where  6h’  is  the  thickness  of  the  radial  catalyst  packed  bed 
and  V/  represents  the  rate  of  production  of  each  species, 
given  in  Eqs.  (18)-(22)  above. 

Eqs.  (13),  (14)  and  (27)  represent  the  complete  coupled 
model  for  fluid  flow  and  reactions  in  the  radial  microreactor, 
which  can  be  solved  using  a  numerical  integration  technique, 
such  as  the  method  of  lines. 

A  similar  analysis  for  the  microchannel  design  (analyzing 
a  single  reactor  microchannel)  gives  the  following  model  for 
the  tubular  configuration,  with  V  representing  the  distance 
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Table  2 

Parameters  used  in  the  simulation — constants 


Cl  =  7.5  x  10  6  molm  2 

Sc  =  102  m2  g  1 

Cl  =  1.5  x  10-5  molm-2 

a 

pb  =  1300  kg  m-3 

Cl  =  7.5  x  10-6  molm-2 

R  =  8.314  J  (mol  K)-1 

Cj  —  1.5  x  10-5  molm-2 

^2a 

T  =  245  °C  =  518. 15K 

Table  3 

Parameters  used  in  the  simulation- 
the  reforming  temperature  of  245  ° 

— k  and  K  parameter  values  calculated  at 
C  =  518  K 

kR  =  31813.7 m2  (mols)-1 
kw  =  86502.8m2  (mols)-1 
kD  =  2732.26m2  (mols)-1 

^H3O(l)=0-68135bar'°'5 

tf*H(|)  =  0.49241  bar-0'5 

^hcoo(I)  =  0-18902  bar 

Kr(  ia)  =  0.59799  bar-1 
Kh(2o)  =  425.41  bar-1 
^ch30(2)  =  3836.6  bar 
^*H(2)  =  3836.6  bar-05 

along  the  reaction  channel  and  tw9  and  ‘f  representing  the 
microchannel  width  and  thickness,  respectively: 


( 


=  m 


( 


m2  m  d M  m  du 

—  — - _l_ - 

P  M  dz  ii  dz 


( 


d 

dz 


=  ( wt)r  j . 


(28) 

(29) 

(30) 


In  order  to  compute  the  reaction  rates  r\  ,  r2  and  r 3,  the 
comprehensive  kinetic  model  developed  by  Peppley  et  al.  [33] 
was  implemented  in  the  simulations.  The  main  equations  of 
this  model  are  given  in  Eqs.  (3 1)— (33)  below  and  the  various 
parameters  used  in  the  model  are  listed  in  Tables  2  and  3. 

All  temperature-dependent  ‘k’  and  4 K 9  values  listed  in 
Table  3  are  calculated  at  the  reforming  temperature  of  245  °C, 
based  on  formulae  available  in  [33]. 


size,  the  radial  design  provides  twice  the  hydrogen  produc¬ 
tion  rate  of  the  microchannel  design,  at  a  fraction  of  the 
pressure  drop.  The  radial  device  fabricated  in  this  work  had 
an  overall  size  (net  reactor  volume)  comparable  to  a  previ¬ 
ously  fabricated  microchannel  based  design  [19]  with  seven 
parallel  microchannels,  each  microchannel  of  dimensions 
400  jutm  x  1000  jam  x  10.4  cm.  Thus,  in  order  to  make  a  fair 
comparison  of  the  expected  pressure  drops,  the  flow  rate  used 
in  the  simulation  of  a  single  microchannel  was  appropriately 
scaled  to  ^th  of  the  total  feed  into  the  microchannel  device. 

As  seen  in  Fig.  8(a),  for  a  methanol  +  water  feed  flow 
rate  sufficient  to  generate  enough  hydrogen  for  a  10-W fuel 
cell ,  the  pressure  drop  across  the  microchannel  device  was 
found  to  be  about  7.8  atm.  In  contrast,  the  pressure  drop 
across  the  radial  device  was  found  to  be  about  0.7  atm  for  a 
methanol  +  water  feed  flow  rate  sufficient  to  generate  enough 
hydrogen  for  a  20 -W fuel  cell  (see  Fig.  9(a)). 

As  seen  in  Fig.  8(b),  the  tubular  (microchannel)  device  has 
a  steep  rise  in  the  flow  velocity  near  the  exit  due  to  the  increase 
in  the  number  of  moles  during  the  reaction  and  expansion  of 
the  reacting  mixture,  which  leads  to  non-uniform  catalyst 
utilization  in  this  region  and  high  pressure  drop  during  the 
operation.  These  effects  are  successfully  countered  in  the  ra¬ 
dial  geometry  and  there  is  a  significantly  lower  variation  in 
the  flow  velocity  magnitude  in  the  radial  device  compared  to 
the  microchannel  device  as  seen  in  Figs.  9(b)  and  8(b),  lead¬ 
ing  to  a  more  uniform  catalyst  utilization  and  significantly 
lower  pressure  drop  in  the  radial  device.  Thus,  the  simulation 
results  support  the  experimental  observations  regarding  the 
higher  throughput  capacity  of  the  radial  device  with  operation 
at  a  fraction  of  the  overall  pressure  drop,  as  noted  in  Section 
5  above. 

It  must  be  noted  that  the  main  objective  of  the  simulations 
was  to  demonstrate  the  advantage  of  the  radial  geometry  in 


kRK CH3CKI)  (/>ch3oh/-Ph25)  0  -  pu2 Pc°2 / kR PCH3OH Ph2o)  C^C^ScPb 

(l  +  Kh  hO(  I )  (^CH3Oh/^h'25)  +  ^7lCOO(  I )  PcOi  +  ^Ol  Id  )  (7h2o/_Ph25))  (l  +  ^ff(ia)^H25) 

kwK*om)(PcoPH2o/P^)(l  -  Pn2Pco2/kwPcoPn2o){CTsfScph 
(1  +  ^CH30(1)  (Ah3Oh/ -Ph25)  +  ^HCOO(l)^C02^H25  +  ^H(1)(^H2o/02 


(31) 


(32) 


fcD^cH3o(2)  (Ah3oh/ ^h25)  (1  -  PlPco/kDPcu3ou)CTS2CTS2aScPb 
(1  +  ^cH30(2)(^>CH30H/f>H25)  +  ^OHap^’O/ 1  +  ) 


where  Cj  ,  ,  Cj  and  Cj  are  the  concentrations  of  ac¬ 

tive  sites  on  the  catalyst  surface,  Sc  is  specific  surface  area 
and  pb  is  catalyst  bulk  density,  respectively  [33]. 

6.2.  Simulation  results 

The  results  from  the  numerical  integration  of  the  above 
models  for  the  radial  and  microchannel  device  are  given  in 
Figs.  8  and  9.  It  is  clear  that  for  the  same  overall  reactor 


providing  a  significant  reduction  in  the  operating  pressure 
drops.  To  that  end,  any  kinetic  model  which  reasonably  rep¬ 
resents  the  reforming  process  would  have  been  adequate.  The 
comprehensive  kinetic  model  [33]  used  in  this  work  is  consid¬ 
ered  to  be  the  most  general  form,  and  was,  therefore,  chosen 
for  the  simulation  studies  to  verify  the  hypotheses  regarding 
the  advantages  of  the  radial  microreactor  configuration. 


A.V.  Pattekar,  M.V.  Kothare  /  Journal  of  Power  Sources  147  (2005)  116-127 


123 


Fig.  8.  (a)  Pressure  profile  in  microchannel  reactor — fluid  pressure  (atm)  vs.  distance  along  the  microchannel  (m);  (b)  velocity  profile  in  microchannel  reactor — 
fluid  velocity  (m/s)  vs.  distance  along  the  microchannel  (m)  (at  a  hydrogen  production  rate  sufficient  for  a  10-W  fuel  cell). 


Fig.  9.  (a)  Pressure  profile  in  radial  reactor — fluid  pressure  (atm)  vs.  distance  from  device  center  (m);  (b)  velocity  profile  in  radial  reactor — fluid  velocity  (m/s) 
vs.  distance  from  device  center  (m)  (at  a  hydrogen  production  rate  sufficient  for  a  20- W  fuel  cell). 


7.  Energy  density  calculations  and  comparisons 

We  will  now  proceed  to  calculate  the  projected  energy 
storage  capacities  possible  using  the  integrated  system  con¬ 
sisting  of  the  radial  microreformer  coupled  with  a  micro  fuel 
cell  and  to  draw  out  important  conclusions  regarding  the  vi¬ 
ability  of  this  technology  in  comparison  with  current  state  of 
the  art  batteries. 

7.7.  The  need  for  alternative  portable  power  sources 

Kelley  et  al.  [34]  and  have  provided  a  comparison  of  the 
specific  energy  density  (cp.  Watt-hour  energy  storage  per  kg 
weight  of  the  battery  system)  possible  using  current  state  of 
the  art  miniature  power  sources  (summarized  in  Table  4). 
Considering  the  projected  improvements  in  the  lithium-ion 
battery  [35],  it  can  be  expected  that  the  energy  storage  ca¬ 
pacity  of  commercially  available  rechargeable  batteries  may 
be  improved  up  to  about  450  Whl-1  or  about  275  Whkg-1 
over  the  next  few  years  [35] .  In  contrast,  based  on  the  DARPA 
palm  power  initiative  guidelines  [36],  it  is  clear  that  several 
key  applications  exist  for  portable  energy  storage  systems 
with  capacities  in  the  range  of  1000  Wh  kg-1  (20  W,  3-h  mis¬ 
sion)  to  over  3000  Whkg-1  (20  W,  10-day  mission).  Thus,  a 
new  micro-power  technology  may  be  the  only  alternative  if 
these  higher  energy  densities  are  to  be  attained  for  portable 
power. 

7.2.  P  EM  fuel  cell  specifications 

For  a  typical  PEM  fuel  cell  operating  at  60%  thermal  ef¬ 
ficiency,  a  hydrogen  flow  rate  of  0.42  g  mol  h-1  would  be 


needed  for  20  W  operation  [37].  Thus,  a  20  h  continuous  op¬ 
eration  at  20  W  would  require  8.40  g  mol  of  H2.  This  trans¬ 
lates  into  a  net  methanol  storage  requirement  of  2.947  g  mol, 
or  94.32  g  methanol,  assuming  2.85  mol  H2  produced  per  mol 
CH3OH  feed  to  the  reformer  as  per  Table  1 ,  for  a  feed  ratio  of 
1.5:1  molar  H20:CH30H.  Now,  a  1.5:1  molar  H20:CH3 OH 
mixture  has  a  density  of  0.9067  g/cm3  at  20  °C  and  contains 
0.49225  g  methanol  per  cc  of  the  mixture  [32].  Thus,  the  net 
water-methanol  (1.5:1  molar)  feed  mixture  storage  require¬ 
ment  for  a  20  W,  20  h  operation  would,  therefore,  be  173.73  g, 
or  191.61  cm3. 

Based  on  the  microfabricated  fuel  cell  sizes  reported 
by  Kelley  et  al.  [34],  a  typical  micro  fuel  cell  for  20  W 
application  (basis:  /zFC  mass  0.031  g,  volume  0.012cm3, 
with  an  active  area  of  0.25  cm2,  operating  at  a  power  density 
of  80  mW  cm-2)  would  have  a  mass  of  0.031  x  20/(0.08  x 
0.25)  =  31  g  and  a  volume  of  12  cm3.  Also,  the  reformer 
reported  in  this  work,  when  scaled  for  the  above  application 
(hydrogen  supply  for  a  20  W  fuel  cell  =  0.42  g  mol  h-1  H2 
production  rate)  would  have  a  mass  of  about  35  g  including 
the  packaging  for  thermal  insulation,  and  a  volume  of  about 
28  cm3 .  Thus,  the  reformer-fuel  cell  combination  would  have 

Table  4 


Specific  energy  density  of  current  state  of  the  art  miniature  power  sources 


Battery  type 

Specific  energy  (Wh  kg  1 ) 

Ni/Cd  battery 

60 

Ni/metal  hydride  battery 

65 

Li/LiAMn204  battery 

130 

LiC6/LixCo02  battery 

90 
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Table  5 


Energy  requirement  for  operation  of  microreformer  producing  hydrogen  for 
a  20  W  fuel  cell 


Reformer  process  step 

Energy  requirement  (W) 

Sensible  heating  of  water  in  feed 

0.6954 

Sensible  heating  of  methanol  in  feed 

0.5504 

Vaporization  requirement  for  water 

2.4483 

Vaporization  requirement  for  methanol 

1.3942 

Endothermic  heat  of  reaction 

2.3850 

Total  power  required  for  operation 

7.4733 

an  intrinsic  mass  of  66  g  and  a  volume  of  about  40  cm3,  in 
addition  to  the  storage  requirements  for  methanol  and  water 
discussed  above.  This  makes  the  total  mass  and  volume  of  the 
power  source  (system  +  fuel)  equal  to  about  (66  +  174)  = 
240  g  and  (40  +  192)  =  232 cm3,  respectively. 

7.3.  Reformer  energy  requirement 

To  make  a  fair  comparison  between  miniature  fuel  cells 
and  state  of  the  art  miniature  batteries,  not  only  does  one 
need  to  consider  a  basis  related  to  the  type  of  application 
being  evaluated  and  the  volume  and  weight  of  the  energy 
conversion  devices  involved,  but  also  the  overall  energy 
requirement  for  operating  the  power  source  (e.g.,  liquid 
feed  vaporization  and  endothermic  reforming  reaction  heat 
requirements,  heat  loss  to  the  ambient,  and  overall  fuel  cell 
efficiency,  etc.  in  the  case  of  the  microreformer-micro  fuel 
cell  combination).  A  detailed  analysis  of  the  the  energy 
requirement  of  each  of  the  process  steps  carried  out  in  the 
microreformer  is  provided  in  Appendix  A  and  summarized 
in  Table  5.  As  shown,  a  reformer  providing  sufficient 
hydrogen  for  a  20  W  fuel  cell  application  would  require 
about  7.4733  W  power  to  vaporize  the  liquid  fuel/water 
mixture  and  carry  out  the  reforming  reaction. 

7.4.  Projected  energy  storage  density 

Assuming  absolutely  no  thermal  integration  between  the 
reformer  and  the  fuel  cell,  if  the  entire  reformer  energy  re¬ 
quirement  is  satisfied  by  electricity  produced  in  the  fuel 
cell,  this  translates  into  a  useful  energy  storage  capacity  of 
20  h  x  (20  -  7.4733)  W  =  250.534  Wh. 

Going  back  to  the  20  W,  20  h  basis,  the  total  mass  and 
volume  of  the  power  source  (fuel  +  system)  are  240  g  and 
232  cm3,  respectively.  This  translates  into  an  energy  storage 
capacity  of  1044  Wh  kg-1  and  1080  Whl-1  for  this  system. 
Of  course,  this  calculation  does  not  consider  the  energy  re¬ 
quirement  for  pumping  the  fuel/water  mixture  at  inlet  pres¬ 
sure,  the  heat  loss  to  the  surrounding  during  the  operation 
and  the  efficiency  of  the  power  conditioning  circuits  that  are 
typically  necessary  in  fuel  cell  based  power  supplies,  which 
would  potentially  bring  down  the  above  values  by  as  much  as 
25 % ,  to  about  783  Wh  kg  - 1  and  about  810  Wh  l- 1 .  However, 
there  is  also  no  allowance  for  heat  integration  between  the  fuel 
cell  and  the  reformer,  which  could  increase  the  overall  sys¬ 


tem  efficiency  and  the  energy  storage  density  by  a  significant 
amount.  For  example,  a  72%  H2  utilization  in  the  fuel  cell 
would  result  in  useful  electrical  energy  output  from  the  fuel 
cell  equal  to  (20  W  x  20  h  x  0.72)  =  288  Wh,  with  a  hydro¬ 
gen  flow  rate  of  0. 1 176  mol  h-1  in  the  fuel  cell  anode  off-gas 
(effluent).  If  the  heating  requirement  for  the  reformer  oper¬ 
ation  is  supplied  by  combustion  of  the  un-utilized  hydrogen 
in  the  anode  off-gas,  this  translates  into  a  heating  capacity 
of  7.8302  W  based  on  the  lower  heating  value  (LHV)  of  hy¬ 
drogen,  which  is  more  than  sufficient  to  provide  the  reformer 
requirement  of  7.4733  W  given  in  Table  5.  This  system  would 
then  have  an  overall  energy  storage  density  of  1 200  Wh  kg- 1 , 
or  about  1241  Wh  1“  1 . 

Recent  developments  in  high  temperature  PEM  fuel  cells 
have  led  to  novel  membranes  that  can  operate  at  high  temper¬ 
atures,  up  to  200  °C  [38].  High  temperature  PEM  fuel  cells 
offer  several  advantages,  including  simpler  water  manage¬ 
ment  and  high  CO  tolerance  (up  to  about  2%  CO  in  the  fuel), 
which  significantly  reduces  the  amount  of  post-reforming 
fuel  clean-up  steps  required.  With  proper  control  of  the  re¬ 
forming  conditions  to  maintain  levels  below  2%,  CO  clean¬ 
up  can  be  completely  avoided,  significantly  simplifying  the 
overall  system.  Apart  from  this,  another  important  advantage 
of  the  high  temperature  PEM  fuel  cell  is  that  the  waste  heat 
is  available  for  exchange  at  a  higher  temperature,  resulting  in 
the  possibility  of  further  heat  integration  in  the  overall  sys¬ 
tem.  For  example,  continuing  the  above  analysis,  40%  of  the 
thermal  energy  in  the  hydrogen  used  in  the  fuel  cell  is  wasted 
as  heat  (60%  thermal  efficiency).  This  translates  into  net  heat 
energy  of  [0.42  m <2/  h~l  x  0.72  (fuel  utilization)  x  0.40  x 
239700.0  J  mol- 1  (LHV  of  H2)]/3600(sh_1)  =  8.054W  at 
200  °C  available  for  further  integration  into  the  system  and 
also  to  take  care  of  the  heat  loss,  providing  a  further  improve¬ 
ment  in  the  energy  storage  density.  Including  a  15%  penalty 
for  the  energy  requirement  related  to  fluid  pumping  and 
power  conditioning  losses,  this  makes  the  overall  energy  stor¬ 
age  densities  of  1200  x  0.85  =  1020  Wh  kg-1,  or  1241  x 
0.85  =  1055  Wh  1_1  easily  attainable  in  a  practical  system. 

Water  produced  in  the  fuel  cell  may  be  recycled  into 
the  system  to  satisfy  the  reformer  requirement,  further 
improving  the  overall  efficiency  and  energy  storage 
capacity  of  the  system.  For  example,  in  a  fuel  cell  op¬ 
erating  at  80%  fuel  utilization,  the  effluent  is  steam  at 
200  °C  at  a  flow  rate  of  (0.42  x  0.80)  =  0.336  g  mol  h_1. 
Noting  that  the  steam  required  in  the  reformer  is 
1.5  x  (0.42/2.85)  =  0.2211  gmolh-1,  the  steam  require¬ 
ment  in  the  reformer  may  be  completely  satisfied  through 
the  fuel  cell  effluent  if  proper  flow  separation  and  control 
systems  are  implemented.  Such  a  highly  integrated  system 
will  also  require  significantly  lower  energy  for  operating  the 
reformer,  since  the  heating  requirements  for  the  water  feed 
are  satisfied  directly  through  the  waste  heat  in  the  fuel  cell. 
For  the  20  W,  20  h  basis,  the  net  volume  and  weight  of  the 
fuel  cartridge  would  then  be  94.32  g,  or  about  121  cm3  (pure 
methanol).  Assuming  a  slightly  larger  fuel  cell  +  reformer 
system  (80  g,  55  cm3)  for  incorporating  the  integration 
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mechanism,  the  net  size  of  fuel  +  system  in  this  case  would 
be  174.32  g  and  176  cm3.  The  net  reformer  energy  require¬ 
ment  would  be  about  4.4268  W  (since  steam  at  200  °C  is 
available  from  the  fuel  cell).  The  hydrogen  flow  in  the  anode 
off-gas  would  be  0.42  x  0.20  =  0.084  mol  h_1,  representing 
a  heating  capacity  of  5.593  W,  which  is  more  than  enough  to 
satisfy  the  reformer  energy  requirement.  The  energy  storage 
capacity  of  this  highly  integrated  system  would  then  be 
(20 W  x  20h  x  0.80)  =  320  Wh.  Again,  including  a  15% 
penalty  for  the  energy  requirement  related  to  fluid  pumping 
and  power  conditioning  losses,  this  translates  into  an  energy 
storage  density  of  1560  Wh  kg-1  or  about  1545  Whl-1,  a 
five-  to  six-fold  improvement  compared  to  the  best  state  of  the 
art  rechargeable  batteries  available  today.  Longer  run  times 
result  in  higher  energy  storage  capacity,  since  the  ratio  of  the 
storage  volume  and  weight  of  the  fuel  to  the  system  would 
increase.  For  a  three-day  run,  this  translates  into  a  stored  en¬ 
ergy  of  (20  W  x  72  h  x  0.80)  =  1152  Wh,  with  a  net  storage 
weight  and  volume  of  (339.55  +  80)  g  and  (435.6  +  55)  cm3, 
respectively.  With  the  15%  penalty  for  system  losses,  this 
would  make  the  energy  storage  density  for  a  3 -day  run  to  be 
2334  Wh  kg-1  or  about  1996  Whl-1.  For  a  10-day  (240  h) 
run,  the  energy  density  would  be  2693  Wh  kg-1  or  about 
2167  Whl-1,  thus,  proving  micro  fuel  cell  technology  to  be 
a  viable  alternative  for  portable  power  needs  of  the  future. 

8.  Conclusion 

A  completely  novel  approach  to  carrying  out  chemical 
reactions  in  a  microreactor  was  presented  in  this  paper.  Con¬ 
ventional  microchannel  based  designs  were  found  to  be  inad¬ 
equate  for  microreaction  systems  which  involve  increase  in 
the  moles  due  to  reaction  and  corresponding  expansion  of  the 
reaction  mixture.  The  radial  device  presented  here  not  only 
solves  the  problems  associated  with  the  microchannel  geom¬ 
etry,  but  also  it  allows  for  rigorous  modeling  and  analysis  of 
the  operation  due  to  the  simple  geometry.  Rigorous  modeling 
and  simulation  studies  of  the  device  operation  revealed  the 
advantages  of  the  new  design  and  confirmed  the  reasoning 
behind  adopting  this  approach  for  developing  the  microreac¬ 
tor.  The  developed  microreactor  produces  sufficient  hydro¬ 
gen  for  fuel  cells  in  the  range  of  up  to  20  W  and  thus  opens 
the  path  for  portable  fuel  cell  based  energy  systems  utilizing 
reforming  technology  for  hydrogen  delivery. 

As  the  demand  for  high  energy  density  power  sources 
grows  with  the  introduction  of  ever  more  power  hungry  mul¬ 
tifunctional  portable  electronic  devices,  miniature  fuel  cells 
are  likely  to  provide  an  excellent  solution  to  the  problems 
faced  by  the  portable  rechargeable  batteries  of  today.  Surely, 
several  challenges  need  to  be  overcome  before  widespread 
introduction  of  this  technology  can  be  successful.  Integra¬ 
tion  of  the  reformer  and  fuel  cell  system,  balance  of  plant 
issues,  efficient  power  conditioning  systems,  and  reliability 
of  long-term  operation  are  just  some  of  the  issues  that  need 
to  be  resolved  in  a  prototype  power  source  utilizing  this  tech¬ 


nology.  As  development  continues  and  these  problems  are 
resolved,  the  microreformer  -  micro  fuel  cell  combination 
with  estimated  energy  storage  densities  in  the  range  of  up 
to  2693  Whkg-1  can  definitely  be  expected  to  appear  in  the 
portable  energy  market  of  the  future. 
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Appendix  A 

A  detailed  analysis  of  the  energy  requirement  of  each  of 
the  process  steps  carried  out  in  the  microreformer  is  provided 
in  this  appendix.  The  results  of  these  calculations  are  sum¬ 
marized  in  Table  5. 

A.l.  Sensible  heating  of  methanol 

The  vapor  pressure  of  methanol  is  given  by  [39]: 
l°gio(^CH3OH) 

3.2447  x  103 

=  45.6171  - - -  13.988  log10(T) 

+  6.6365  x  10“3r  -  1.0507  x  10"13T2  (A.l) 

where  Pq+oh  *s  measured  in  mmHg  and  T  is  the  absolute 
temperature.  Thus,  at  the  reformer  inlet  pressure  of  1.7  atm 
(see  Fig.  9(a)),  the  boiling  point  of  CH3OH  is  351.85  K,  ob¬ 
tained  from  the  above  equation. 

The  specific  heat  capacity  of  liquid  methanol  C™3°H,(P 

(J  (mol  K)_1)  is  given  by  [39]: 

cch3oh,(D  =  40.152  +  3.1046  x  10_1T  -  1.0291 

x  10"3r2  +  1.4598  x  10"6r3  (A. 2) 

where  T  is  the  absolute  temperature.  Therefore,  the  sensible 
heat  requirement  for  liquid  methanol  from  a  room  tempera¬ 
ture  of  298. 15  K  to  its  boiling  point  of  35 1 .85  K  is: 

>-351.85 

/  C+OH  (1)d7’  =  4433.458  J  mol-1.  (A.3) 

7298.15 

Also,  the  specific  heat  capacity  of  methanol  vapor  C™30H,(^ 
(J  (mol  K)_1)  is  given  by  [39]: 
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cCH3OH,(g)  =  40>046  _  3.8287  x  10 ~2T  +  2.4529 

x  10_4r2  -  2.1679  x  10-77’3  +  5.99 
x  10-117’4.  (A.4) 


Therefore,  the  sensible  heat  requirement  for  methanol  vapor 
from  its  boiling  point  of  351.85  K  to  the  reactor  operating 
temperature  of  245  °C  =  518.15  K  is: 


/*518.15 
J  351.85 


cCH3OH,(g)dr  =  9009.682  J  mol-1. 


Thus,  the  total  sensible  heat  load  for  methanol  is: 
AR+3°H’sensible  =  13443. 14  J  mol-1. 


(A. 5) 


(A. 6) 


For  a  methanol  flow  rate  of  0. 1474  mol  h- 1  (sufficient  to  pro¬ 
duce  0.1474  x  2.85  =  0.420  mol  h-1  of  H2),  this  translates 
into  an  energy  requirement  of  13443.14  x  0.1474/3600  = 
0.5504  W. 


A.  2.  Methanol  vaporization 


The  enthalpy  of  vaporization  of  methanol  is  given  by  [39] : 


A//ch3oh 


52.723 


V  512.58 


0.377 

kJ  mol-1. 


(A. 7) 


Thus,  the  methanol  vaporization  heating  requirement  at  the 
boiling  temperature  of  351.85  K  is 


52.723 


(  351.85 

V  ~~  512.58 


0.377 

kJmol-1  x 


0.1474 

3600 


mol  s  1 


=  1.3942  W. 


(A. 8) 


A3.  Sensible  heating  of  water 


Therefore,  the  sensible  heat  requirement  for  steam  from  its 
boiling  point  of  388.70  K  to  the  reactor  operating  temperature 
of  245  °C  =  518.15  K  is: 

/*518.15 

/  c“2O(g)dT  =  4501.367Jmol-1.  (A.12) 

7 388.70 

Thus,  the  total  sensible  heat  load  for  water  is: 

AHtota?’SenSible  =  11322.71  Jmol-1.  (A.13) 

For  a  water  flow  rate  of  0.22 1 1  mol  h- 1  (sufficient  to  produce 
(0.2211/1.5)  x  2.85  =  0.420 molh-1  of  H2),  this  translates 
into  an  energy  requirement  of  11322.71  x  0.2211/3600  = 
0.6954  W. 

A.4.  Water  vaporization 

The  specific  enthalpy  of  saturated  water  (liquid)  at 

388.70  K  is  8727.15  J  mol-1,  and  that  of  saturated  steam  (va¬ 
por)  is  48591.65  J  mol-1  [32].  Thus,  the  energy  requirement 
for  vaporization  of  the  feed  water  at  the  microreformer  inlet 
pressure  of  1.7  atm  is 

,  0.2211  , 

(48591.65  -  8727. 15) Jmol-1  x  - mols-1 

3600 

=  2.4483  W.  (A.  14) 


A.  5.  Endothermic  heat  of  reaction 

Considering  the  overall  reforming  reaction 

CH3OH(g)  +  H20(g)  =  CO{f  +  3H71  (A.  15) 

we  have: 


At  the  reactor  feed  pressure  of  1.7  atm,  the  boiling  point 
of  water  is  388.70  K,  obtained  from  steam  tables  [32].  For 
calculating  the  liquid  water  heating  requirement,  we  have  the 
following  correlation  for  specific  heat  capacity  of  water  [39]: 

ch2o,(1)  =  92.053  -  3.9953  x  10-2T  -  2.1103 

x  10-4T2  +  5.3469  x  10-7r3  J(molK)-1(A.9) 

Therefore,  the  sensible  heat  requirement  for  liquid  water 
from  a  room  temperature  of  298.15  K  to  its  boiling  point 
of  388.70  K  is: 

/*388.70 

/  Cp2°,(1)d7’  =  6821.339  Jmol-1.  (A.10) 

7298.15 


AT/298.15  at  =  48.97  kJmol  1 


(A.  16) 


>  A //245  °C  =  A //518J5  k 


=  A//298.15  K  + 


/*  518.15 
7298.15 


ACp(T)dT 


=  58. 25  kJmol  1 . 


(A.  17) 


where 


^  c  pco2 


rte)  -  r(g) 

Ph2o  /?ch3oh 


(A.  18) 


which  is  given  by  [39]: 


Also,  the  specific  heat  capacity  of  water  vapor  C^2°,(^ 
(J  (mol  K)-1)  is  given  by  [39]: 

Cjk°’(g)  =  33.933  _  8.4186  x  10-3T  +  2.9906 

x  10-5T2  -  1.7825  x  io-8r3 
+  3.6934  x  10-12T4.  (A.l  1) 


A Cp  =  29.655  +  0.149555T  -  4.10398  x  10-4T2 
+  3.34252  x  10-7T3  -  9.01766  x  10-11T4 

(A.  19) 

Thus,  the  heat  requirement  for  carrying  out  the  endothermic 
reforming  reaction  is  58250  x  0.1474/3600  =  2.3850  W. 
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Appendix  B.  Supplementary  data 

Supplementary  data  associated  with  this  article  can 
be  found,  in  the  online  version,  at  10.1016/j.jpowsour. 
2005.01.024. 


References 

[1]  P.  McFadden,  BIOSENSORS:  broadband  biodetection:  holmes  on  a 
chip,  Science  297  (2002)  2075. 

[2]  V.  Percec,  Chemistry:  catalyst  hunt  accelerates,  Nature  424  (2003)  135. 

[3]  T.R.  Boussie,  G.M.  Diamond,  C.  Goh,  K.A.  Hall,  A.M.  LaPointe,  M. 
Leclerc,  C.  Lund,  V.  Murphy,  J.A.W.  Shoemaker,  U.  Tracht,  H.  Turner, 
J.  Zhang,  T.  Uno,  R.K.  Rosen,  J.C.  Stevens,  A  fully  integrated  high- 
throughput  screening  methodology  for  the  discovery  of  new  polyolefin 
catalysts:  discovery  of  a  new  class  of  high  temperature  single-site  group 
(IV)  copolymerization  catalysts,  J.  Am.  Chem.  Soc.  125  (2003)  4306- 
4317. 

[4]  M.A.  Burns,  B.N.  Johnson,  S.N.  Brahmasandra,  K.  Handique,  J.R. 
Webster,  M.  Krishnan,  T.S.  Sammarco,  PM.  Man,  D.  Jones,  D. 
Heldsinger,  C.H.  Mastrangelo,  D.T.  Burke,  An  integrated  nanoliter 
DNA  analysis  device,  Science  282  (1998)  484-487. 

[5]  R.F.  Service,  Shrinking  fuel  cells  promise  power  in  your  pocket,  Sci¬ 
ence  296  (2002)  1222. 

[6]  S.  Park,  J.M.  Vohs,  R.J.  Gorte,  Direct  oxidation  of  hydrocarbons  in  a 
solid-oxide  fuel  cell,  Nature  404  (2000)  265-267. 

[7]  G.A.  Deluga,  J.R.  Salge,  L.D.  Schmidt,  X.E.  Verykios,  Renewable  hy¬ 
drogen  from  ethanol  by  autothermal  reforming,  Science  303  (2004) 
993-997. 

[8]  L.  Schlapbach,  A.  Zuttel,  Hydrogen-storage  materials  for  mobile  ap¬ 
plications,  Nature  414  (2001)  353. 

[9]  U.  Biinger,  W.  Zittel,  Hydrogen  storage  in  carbon  nanostructures — 
still  a  long  road  from  science  to  commerce?  Appl.  Phys.  A  72  (2001) 
147-151. 

[10]  P.  Pfeifer,  M.  Fichtner,  K.  Schubert,  M.A.  Liauw,  G.  Emig,  Mi- 
crostructured  catalysts  for  methanol-steam  reforming,  in:  Proceed¬ 
ings  of  the  3rd  International  Conference  on  Microreaction  Technology, 
DECHEMA,  Frankfurt,  Germany,  1999. 

[1 1]  A.Y.  Tonkovich,  J.L.  Zilka,  M.J.  LaMont,  Y.  Wang,  R.S.  Wegeng,  Mi¬ 
croChannel  reactors  for  fuel  processing  applications.  I.  Water  gas  shift 
reaction,  Chem.  Eng.  Sci.  54  (1999)  2947-2951. 

[12]  P.J.  de  Wild,  M.J.F.M.  Verhaak,  Catalytic  production  of  hydrogen  from 
methanol,  Catal.  Today  60  1  (2000)  3-10. 

[13]  S.  Fitzgerald,  R.  Wegeng,  A.  Tonkovich,  Y.  Wang,  H.  Freeman,  J. 
Marco,  G.  Roberts,  D.  VanderWeil,  A  compact  steam  reforming  re¬ 
actor  for  use  in  an  automotive  fuel  processor,  in:  R.S.W.W.  Ehrfeld, 
I.  Rinard  (Eds.),  Proceedings  of  the  4th  International  Conference  on 
Microreaction  Technology,  Atlanta,  GA,  2000,  pp.  358-363. 

[14]  A.V.  Pattekar,  M.V.  Kothare,  S.V.  Karnik,  M.K.  Hatalis,  A  microreac¬ 
tor  for  in  situ  hydrogen  production  by  catalytic  methanol  reforming, 
in:  Proceedings  of  the  5th  International  Conference  on  Microreaction 
Technology  (IMRET  5),  Strasbourg,  France,  2001. 

[15]  J.D.  Holladay,  E.O.  Jones,  M.  Phelps,  J.  Hu,  Microfuel  processor  for 
use  in  a  miniature  power  supply,  J.  Power  Sources  108  (2002)  21. 

[16]  A.V.  Pattekar,  M.V.  Kothare,  Design  and  fabrication  of  a  microreactor 
for  hydrogen  production  by  catalytic  methanol  reforming,  in:  Proceed¬ 
ings  of  the  2002  AIChE  Annual  Meeting,  Indianapolis,  IN,  2002. 

[17]  M.V.  Kothare,  A.V.  Pattekar,  Fuel  processing  microreactors  for  hydro¬ 
gen  production  by  methanol  reforming,  in:  Proceedings  of  the  Knowl¬ 
edge  Foundation’s  5th  International  Small  Fuel  Cells  2003:  Small  Fuel 
Cells  for  Portable  Power  Applications,  New  Orleans,  LA,  2003. 


[18]  A.V.  Pattekar,  M.V.  Kothare,  Microreactor  and  method  of  use  to  pro¬ 
duce  hydrogen  by  methanol  reforming,  US  Patent  Pub.  App.  no. 
20040179980.  Provisional  disclosure  filed  December,  2002  (in  final 
form:  December  2003). 

[19]  A.V.  Pattekar,  M.V.  Kothare,  A  microreactor  for  hydrogen  production 
in  micro  fuel  cell  applications,  J.  Microelectromech.  Syst.  13  (2004) 
7-18. 

[20]  M.  Richter,  A.  Drost,  J.  Kruckow,  A  high  performance  silicon  microp¬ 
ump  for  fuel  handling  in  DMFC  systems,  in:  Proceedings  of  the  2003 
Fuel  Cell  Seminar,  Miami  beach,  FL,  2003. 

[21]  M.A.  Unger,  H.  Chou,  T.  Thorsen,  A.  Scherer,  S.R.  Quake,  Monolithic 
microfabricated  valves  and  pumps  by  multilayer  soft  lithography,  Sci¬ 
ence  288  (2000)  113-116. 

[22]  W.  Qu,  I.  Mudawar,  Measurement  and  prediction  of  pressure  drop  in 
two-phase  micro-channel  heat  sinks,  Int.  J.  Heat  Mass  Transfer  46 
(2003)  2737-2753. 

[23]  A.A.  Savoretti,  D.O.  Borio,  V.  Bucala,  J.A.  Porras,  Non-adiabatic 
radial-flow  reactor  for  styrene  production,  Chem.  Eng.  Sci.  54  (1999) 
205-213. 

[24]  A.V.  Pattekar,  M.V.  Kothare,  Novel  microfluidic  interconnectors  for 
high  temperature  and  pressure  applications,  J.  Micromech.  Microeng. 
13  (2003)  337-345. 

[25]  See  http://www.lehigh.edu/~asp2/mov/movie-l-uni-flow.wmv  for  a 
video  of  the  uniform  flow  distribution  of  the  fluid  exiting  the  vaporiza¬ 
tion  region  (also  available  as  part  of  supplementary  online  material  at 
the  Journal  of  Power  Sources  web-site). 

[26]  See  http://www.lehigh.edu/~asp2/mov/150C-various-flow-rates.wmv 
for  a  video  of  vaporizer  operation  at  150  °C  and  various  liquid  water 
feed  flow  rates.  (Also  available  as  part  of  supplementary  online  material 
at  the  Journal  of  Power  Sources  web-site). 

[27]  See  http://www.lehigh.edu/~asp2/mov/180C-various-flow-rates.wmv 
for  a  video  of  vaporizer  operation  at  180  °C  and  various  liquid  water 
feed  flow  rates.  (Also  available  as  part  of  supplementary  online  material 
at  the  Journal  of  Power  Sources  web-site). 

[28]  See  http://www.lehigh.edu/~asp2/mov/200C-various-flow-rates.wmv 
for  a  video  of  vaporizer  operation  at  200  °C  and  various  liquid  water 
feed  flow  rates.  (Also  available  as  part  of  supplementary  online  material 
at  the  Journal  of  Power  Sources  web-site). 

[29]  J.C.  Amphlett,  M.J.  Evans,  R.A.  Jones,  R.F.  Mann,  R.D.  Weir,  Hydro¬ 
gen  production  by  the  catalytic  steam  reforming  of  methanol.  Part  1. 
The  thermodynamics,  Can.  J.  Chem.  Eng.  59  (1981)  720-727. 

[30]  W.L.  McCabe,  J.C.  Smith,  P.  Harriot,  Unit  Operations  of  Chemical 
Engineering,  fifth  ed.,  McGraw-Hill,  New  York,  1993  (Chapter  7). 

[31]  R.B.  Bird,  W.E.  Stewart,  E.N.  Lightfoot,  Transport  Phenomena,  Wiley, 
New  York,  NY,  1960. 

[32]  R.  Perry,  D.  Green,  Chemical  Engineer’s  Handbook,  sixth  ed., 
McGraw-Hill,  New  York,  NY,  1984. 

[33]  B.A.  Peppley,  J.C.  Amphlett,  L.M.  Kearns,  R.F.  Mann,  Methanol  steam 
reforming  on  Cu/Zn0/Al203  catalysts.  Part  2.  A  comprehensive  kinetic 
model,  Appl.  Catal.  A:  Gen.  179  (1999)  31-49. 

[34]  S.  Kelley,  G.  Deluga,  W.  Smyrl,  Miniature  fuel  cells  fabricated  on 
silicon  substrates,  AIChE  J.  48  (2002)  1071-1082. 

[35]  A.V.  Pattekar,  M.V.  Kothare,  Miniature  fuel  processors  and  fuel  cells  for 
portable  power:  industry  trends  and  potential  applications,  in:  Proceed¬ 
ings  of  the  2003  AIChE  Annual  Meeting,  San  Francisco,  CA,  2003. 

[36]  See  http://www.darpa.mil/dso/thrust/matdev/palmpower/research.html 
for  a  summary  of  the  goals  of  the  DARPA  palm  power  project. 

[37]  Fuel  Cell  Handbook,  sixth  ed.,  EG&G  Services  Inc.,  US  Department 
of  Energy,  Morgantown,  WV,  November  2002. 

[38]  A.  Schechter,  R.F.  Savinell,  Imidazole  and  1 -methyl  imidazole  in  phos¬ 
phoric  acid  doped  polybenzimidazole,  electrolyte  for  fuel  cells,  Solid 
State  Ionics  147  (2002)  181-187. 

[39]  C.L.  Yaws  (Ed.),  Chemical  Properties  Handbook,  McGraw-Hill,  New 
York,  NY,  1999. 


